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Rokhsar-Kivelson Models of Bosonic Symmetry-Protected Topological States 
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Perimeter Institute for Theoretical Physics, Waterloo, ON, N2L 2Y5, Canada 

A platform for constructing microscopic Hamiltonians describing bosonic symmetry-protected topological 
(SPT) states is presented. The Hamiltonians we consider are examples of frustration-free Rokhsar-Kivelson 
models, which are known to be in one-to-one correspondence with classical stochastic systems in the same 
spatial dimensionality. By exploring this classical-quantum mapping, we are able to construct a large class 
of microscopic models which, in a closed manifold, have a non-degenerate gapped symmetric ground state 
describing the universal properties of SPT states. Examples of one and two dimensional SPT states which 
illustrate our approach are discussed. 


I. INTRODUCTION 


The SO (3) symmetric spin-1 anti-ferromagnetic Heisen¬ 
berg chain, which was shown by Haldane 1 3 to have a symme¬ 
try preserving gapped ground state, provides the oldest known 
example of a bosonic symmetry-protected topological (SPT) 
state in one dimension. The topological character of this state 
is captured by a topological 0-term present in the non-linear 
sigma model effective action describing long-wavelength de¬ 
grees of freedom. 1 3 

An insightful account of the properties of this anti¬ 
ferromagnetic spin chain was given by Affleck, Kennedy, 
Lieb, and Tasaki (AKLT), who constructed a Hamiltonian in 
the same phase as the Heisenberg model, where the S = 1 
spins emerge from the composition of underlying 5=1/2 
degrees of freedom.® 5 With periodic boundary conditions, the 
AKLT model has a non-degenerate ground state that does not 
break any symmetries and is separated from the first excited 
state by a finite gap. With open boundary conditions, on the 
other hand, the AKLT model makes it manifest that each edge 
supports a “free” 5=1/2 degree of freedom contributing to 
a 2-fold degeneracy per edge. Interestingly, while the Hamil¬ 
tonian is 50(3) symmetric and the bulk degrees of freedom 
transform linearly under this symmetry, the effective 5 = 1/2 
spins on the edges transform projectively under the action of 
the spin rotation symmetry: an initial 5 = 1/2 spin state ro¬ 
tated by 360 degrees about an arbitrary axis is mapped into 
itself up to a minus sign. The inter-connection among the 
topological 0-term action, the ground state degeneracy with 
open boundary conditions and the projective representation of 
the global symmetry on the edge degrees of freedom makes 
this system a non-trivial gapped phase of matter. 

Inspired by the example of the anti-ferromagnetic chain, 
there have been recent proposals to classify gapped SPT 
phases of matter protected by a global symmetry G using var¬ 
ious mathematical frameworks such as group cohomology® 
which generalizes the concept of projective representations, 
topological field theories^*® and non-linear sigma models in 
the presence of a topological 0-term action compatible with 
the global symmetry G. 1718 Recently, a number of micro¬ 
scopic models of bosonic SPT states have been studied, which 
help to shed light on the role played by physical interactions 
in bringing about SPT phases PBSalsl Mg lt^ 

The purpose of this paper is to provide a framework for con¬ 
structing microscopic models capable of describing bosonic 


SPT states. As we shall see, some of the exactly solvable 
models previously studied 2021 ® will be identified as special 
cases of a large class of models to be constructed here. We 
shall also be able to construct parent Hamiltonians for two di¬ 
mensional Z 2 x Z 2 paramagnets, whose effective edge theory 
was shown in Ref. [16| to be in direct relation to non-trivial 
3-cocycles. 

The classes of gapped bosonic insulators protected a by 
global symmetry G that we shall be concerned with have, 
on a d-dimensional closed manifold, a non-degenerate ground 
state 


I ) = v7 2= E \s), ( 1 . 1 ) 

where {| s )} denotes an orthonormal many-body basis, 
Eg{s) € K. is a non-universal local function related to the 
decay of correlations of local operators in the ground state 
and the phase Wg(s) € K. is a universal piece that endows the 
ground state O* with its non-trivial topological properties. 
Wg(s) plays, at the microscopic level considered here, a role 
analogous to the topological 0-term.® 1718 

When Wg{s) = 0, one obtains from Eq. CD the nodeless 
ground state 


I $g) = 


Vm 


E 


-§ E G (s) 


s). 


( 1 . 2 ) 


The form of the ground state Eq. (1.2 1 is very appealing; for 
equal time correlation functions of operators in the diagonal 
representation {| s )}, 


; $ G | d a (§) 6 b (§) | $ G ) = E 


„-ftE G (s) 

z(l3) 


(1.3) 


can be interpreted as equal time correlation functions of an 
equilibrium d-dimensional statistical mechanical system with 
classical configurations {s}, each one occurring with prob¬ 
ability p i 0 ' 1 = e~^ E g ^/Z{ j5), where the real parameter /3 
acquires the natural interpretation of an effective inverse tem¬ 
perature and the normalization factor of the ground state, 

Z(« = ]>>-»*«<•>, „. 4) 
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is interpreted as the partition function of the classical system 
in “thermal” equilibrium. Hence, if the associated classical 
model described by the partition function Eq. ( |1.4| > is in the 
“disordered” phase, then typical correlation functions of lo¬ 
cal operators distant by |r| behave as l^l/■£, for some finite 
correlation length £, and the representation (1.21 can be asso¬ 
ciated with a quantum many-body ground state in its gapped 
phase. 

In fact the foregoing classical -qua ntum correspondence can 
be made more precise. 27 31 Eq. (1.21 is recognized as the zero 
energy ground state of a class of quantum dimer-like models 
at the so-called Rokhsar-Kivelson (RK) point. 27 In Ref. [27] it 
was noted that dimer-dimer correlation functions of the square 
lattice quantum dimer model at the RK point can be computed 
exactly from the corresponding classical dimer problem. 32 In 
Ref. |29j Ardonne, Fendley and Fradkin have established the 
quantum-classical correspondence to more general classes of 
RK Hamiltonians beyond dimer models. In Ref. [30] Henley 
has observed that any stochastic classical system described 
by a real transition rate matrix M can be interpreted, via a 
similarity transformation, as an RK Hamiltonian. In Ref. [3T1 
Castelnovo, Chamon, Mudry and Pujol have shown that the 
reverse is true, namely, that given a quantum RK Hamiltonian 
in a “preferred basis” [the basis in which the ground state is 
expressed as a linear combination of same-sign coefficients 
as in Eq. ( | 1 .2[ >], there exists an associated stochastic classical 
model whose spectrum of relaxation rates is (up to an overall 
minus sign) the same as the energy spectrum of the quantum 
RK Hamiltonian and whose equilibrium probability distribu¬ 
tion is the square of the coefficients in the expansion of the 
RK quantum ground state Eq. 0 - 

In light of the above arguments, if one considers the con¬ 
figurations {s} to be made of spins on a lattice, then the 
ground state Eq. 0 offers a natural representation of a para¬ 
magnetic state, provided the corresponding classical system is 
chosen to have a spectrum of relaxation rates with a finite gap 
and correlation functions of local operators, Eq. ( |1.3[ ), exhibit¬ 
ing short-range behavior. 

As for the role played by symmetries, we now let the quan¬ 
tum system be invariant under a global symmetry group G, 
whose action on the basis | s ) is represented by 


hence. 


S G \s) = \gs). 


(1-5) 


It is then clear that the ground state Eq. 0 is a unique and 
G invariant state provided the local “classical energy” E G (s) 
is symmetric under the transformation {s} —► {5 s}: 


^g) =W g |$ g ) 


(1.7b) 


For the SPT ground state Eq. 0 to be invariant under the 
symmetry G, it is required that 


Eg{9s) = E g (s) , 

W G (g s ) = W G {s) mod 2-7T , 


(1.8a) 


if G is an unitary global symmetry, and 


E G {gs) = E g {s) , 

W G (g s ) = — W G {s ) mod 27 t , 


(1.8b) 


if G is an anti-unitary global symmetry. 

Now let H g and H G be, respectively, the quantum Hamil¬ 
tonians whose non-degenerate ground states are | ( l> ( ;) and 
| 'Ifo;). Then the unitary mapping Eq. (1.7b 1 establishes 


H g = W g Hg W q 1 


(1-9) 


Starting, thus, from a parent Hamiltonian H G for the trivial 
gapped state Eq. s one can construct a parent Hamilto¬ 
nian H G for the SPT ground state Eq. 0 using Eq. ( |1.9[ ), if 
the unitary transformation connecting the two ground states, 
Eq. ( |1.7[ i, is known. 

This paper is organized as follows. In Sec. 0 we review 
the relevant points about the mapping between stochastic clas¬ 


sical systems and quantum RK Hamiltonians. In Sec. (III 1 


we construct parent Hamiltonians of one dimensional bosonic 
SPT states with Z n x Z„ symmetry, where we shall use 
the concept of entanglement spectrum degeneracy to deter¬ 
miner the unitary transformation Eq. 0 that maps the trivial 
ground state to other n — 1 topological phases. In Sec. ( |IV| ) we 
discuss the one dimensional SPT state with anti-unitary time- 
reversal symmetry Z]f. In Sec. ([vji we construct two dimen¬ 
sional microscopic models that account for all the 8 possible 
classes of Z2 XZ2 SPT paramagnets. We show that the Z2 XZ2 
symmetry transformation projected onto the one dimensional 
edge acquires a non-onsite form, which was studied in Ref. [16] 
in connection with non-trivial 3-cocycles in the group coho¬ 
mology. Finally we draw some conclusions and point to future 
directions in Sec. m- 


E G (gs ) = E G (s). 


( 1 . 6 ) 


II. CLASSICAL-QUANTUM MAPPING 


The central point of this paper is the observation, which 
will be supported by concrete examples, that in a closed man¬ 
ifold, the SPT ground state Eq. ( |1.1[> c an be obtained from the 
trivial insulator ground state Eq. ( |1.2) via a global symmetry¬ 
preserving unitary transformation W G , whose action on the 
many-body basis {| s )} is 

W G |s) :=e iWc(s) | a), (1.7a) 


In order to make the discussion self-contained, we review 
the essential aspects of the relationship between quantum RK 
Hamiltonians and stochastic classical systems.®® The con¬ 
tent of Secs. ( |II A[ > and ( |IIB| > closely follows Ref. [3j which 
the reader may consult for further details. In Sec. ( jlTC) we 
give the general form of the SPT-RK Hamiltonians describing 
bosonic SPT states in d-dimensional space obtained from the 
unitary mapping Eq. (1.91. 
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A. Rokhsar-Kivelson Hamiltonians 


A quantum RK Hamiltonian satisfies three propertied- 31 : 

1. The orthonormal elements of the basis E = {| s)}, 
which span the Hilbert space, 

(s|s') = 6 S , S >, ll=^|s)(s|, (2.1) 

S 

form a countable set. 


2. The quantum Hamiltonian can be decomposed into a 
sum of positive-semidefinite projector-like Hermitian 
operators P s s i [Eq. (2.21], 


3. The ground state [Eq. (2.41] is annihilated by every 
P s s ' and the normalization constant of the ground state 
can be interpreted as the partition function of a classical 
system [Eq. 


The RK Hamiltonian takes the form 

PrK = .. 'y ^s,s f Ps,s' 7 
(«>«') 

where 


^ ^ 7 ^s,s' P 0 7 


(2.2a) 


(2.2b) 


Ps,s'=~\s)(s'\-\S’)(S\ 

+ e-« [*<«')-*(«)] | s )< s | + e-f | a > )( s ' |, 

(2.2c) 

where j3 £ R, P(s) £ K and P S)S / is a projector-like positive 
semi-definite Hermitian operator satisfying 

P S 2 S , = 2 cosh | ^ [P(s') - E(s)] | P 8 ,y ■ (2.3) 


One easily verifies that the nodeless state 
1 


with normalization constant 


E 


e -f 


m = E 


s) 7 


-fiE(s) 


(2.4) 


(2.5) 


satisfies 


P, 


‘J’RK ) — 0 , 


V (s,s'). 


( 2 . 6 ) 


Thus E q. \2A) is the ground state of the RK Hamilto¬ 
nian Eq. (2.21 with energy zero. The normalization constant 


Eq. ( |2.5| > can be interpreted as the partition function of a clas¬ 
sical system with classical energy E(s) at the effective inverse 
temperature 0. In the “infinite temperature” limit 0 = 0, the 
Boltzmann factors in the partition function tend to unity and 
Z( 0) counts the number of allowed configurations. 


B. Relation between quantum RK Hamiltonians and 
stochastic classical systems 

In order to show the relation between quantum RK Hamil¬ 
tonian and stochastic classical systems, one considers the real¬ 
valued matrix M defined by 30 31 

Mss' = -e~l (P RK ) SS , , (2.7) 

where (H RR ) SS , = (s | Prk | s ') denotes the matrix elements 
of I Irk- From the fact that I Irk is a real and Hermitian oper¬ 
ator, it follows that 


(Hrk) ss , = (H rk ) s , s (2.8) 

and, from Eqs. \2.2) and ( |2.7[ ), that the matrix M satisfies: 

M ss ' > 0 , if sjts', (2.9a) 


and 


Ms, 


s'^S 


Y ■ 


(2.9b) 


One then considers a classical system with a phase space 
formed by configurations {s}, which, as a function of time r, 
can be visited stochastically with probability p s (r) evolving 
according to the master equation 


dpsjr) 

dr 


= Y M ss'Ps'{t) 


s' 


E \ M ss'Ps'(t) - M s 'sPs{t) 


( 2 . 10 ) 


where Eq. < 2.9b) has been used to achieve the last equality 
of Eq. (2.10i. The first term on the r.h.s. accounts for the 
transitions out of the configurations s' into the configuration 
s, while the second term on the r.h.s. accounts for transitions 
out of the configuration s into the configurations s'. 

Moreover Eqs. ( |2.7| and ( |2.8[ ) are easily seen to imply, for 
every pair of indices (s, s'), 


Mss' P^ = Ms'sP, 


( 0 ) 


(2.11a) 


where 


p(°) = — e 

Ps ~m 


-p E(s) 


(2.11b) 


Eq. (2.11 1> implies the condition of detailed balance on the 


matrix M as well as that pi 0 -* is the equilibrium probability 


distribution associated with the classical dynamics Eq. (2.10|>. 


Denoting by X n and t/>, 




respectively, the right- 


eigenvalues and right-eigenvectors of M, 


Y M ss' = A„ ipi R ' n) , (2.12) 
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then the time dependent solution of Eq. (2.10) can be ex¬ 
pressed as 

Ps(t) = a n (0)e XnT 'tp { s R ’ n) , (2.13) 


largest characteristic time scale associated with the decay into 
the equilibrium configuration is finite in the thermodynamic 


limit, then it follows from Eq. (2.14) that the many-body en¬ 


ergy spectrum of the quantum Hamiltonian possess a finite 
energy gap for excitations above ground state. 


where a n (0) are coefficients determined by the initial condi¬ 
tions. 

Since Eq. ( |2.7j i establishes, up to an overall minus sign, a 
similarity transformation between // KK and M, the spectrum 
of relaxation rates {A„ } of M and the energy spectrum {e n } 
of Ark are simply related: 

e n = —X n > 0. (2.14) 


When the classical system whose dynamics is described by 
Eq. (|2.10[) has a spectrum of relaxation rates such that the 


C. SPT-RK Hamiltonians 


We now give the general form of the SPT-RK Hamiltonians 
with ground state given by Eq. CD- 

Let the d-dimensional quantum system, protected by global 
symmetry G, in its trivial phase be described by an RK Hamil¬ 
tonian Hq of the form Eq. \2.2\ with the non-degenerate 
ground state | $g), Eq. Q where we impose the symmetry 
constraint Eq. (1.6) upon Eq{s). Then the unitary transfor¬ 
mation Eq. (1.71 yields the SPT ground state | )> Eq. CD, 

and Eq. ( |1.9| > yields the SPT-RK Hamiltonian He'. 


HLg = W g Hg W G ' 


1 

2 


s'^S 

w s>s' tv > 

(s,s') 


(2.15a) 


where 


w s ,yeR, lo s , s '> 0, (2.15b) 

Vs,s' = — e* [WeG)-W G (s')] | g g' | _ e i [Wo.(«')-Wo(«)] | s' ^ g | + e -| [Ba^-BoW] | 8 ^ 8 | + e ~§ [® G («)-Eo(«')] | s ' )( s ' | . 

(2.15c) 


Before we proceed to discuss specific SPT systems, we 
close this section with a few important remarks: 

(a) In Secs. ( |H A| ) and ( |ll B[ i we saw that the structure of the 
RK ground state is determined only by ft E(s), while 
the energy spect rum is determined solely by the cou¬ 
plings w SjS 'p2ED 


(b) If the couplings are such that the energy spec¬ 
trum has a non-degenerate gapped ground state (which 
is the only case considered in this work), then Eq. ( |2.4| > 
parametrizes a large class of ground states, whereby 
P E(s) controls the correlation length associated with 
the decay of correlation functions such as Eq. (1.31. 


(c) In the particular case ft = 0, the ground state Eq. ( |2.4[ > 
acquires the form of a trivial product state. It will 
prove useful to work in this “infinite temperature” limit, 
where the correlation length is zero, in order to deter¬ 
mine the form of the unitary transformation Eq. CD- 


(d) While the unitary mapping Eq. CD leaves the energy 
spectrum unchanged, it introduces a non-trivial entan¬ 
glement structure in the SPT ground state Eq. CD that 
is responsible for its universal topological properties. 


III. SPT STATES IN ONE DIMENSION 


In this section we derive parent Hamiltonians of one dimen¬ 
sional bosonic SPT states protected by Z n x Z„ symmetry. 
In a chain with periodic boundary conditions, each of these 
phases is described by a non-degenerate gapped symmetric 
ground state. In a chain with open boundary conditions, on the 
other hand, there remains a trivial phase with a non-degenerate 
ground state, while n — 1 phases have n-fold degeneracy per 
edge, accounting for a total n 2 -fold degeneracy of the ground 
state manifold in the thermodynamic limit. 


Our aim is to construct the unitary transformations Eq. ( |1.7[ ) 
connecting the trivial ground state and the n — 1 non-trivial 
SPT ground states. Our strategy in deriving such unitary map¬ 
pings is to draw on the notion of entanglement spectrum de¬ 
generacy as follows. As we pointed out in the remark (c) 
of Sec. (II Ci, in the “infinite temperature” limit /? = 0 the 
ground state of the trivial SPT chain reduces to a product 
state. The entanglement structure of the product state is as 
simple as it gets, for the Schmidt decomposition with respect 
to any partition contains only one eigenvalue (equal to 1). We 
shall find the unitary transformation Eq. CD by demanding 
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that, in the limit (3 = 0, the entanglement spectrum of the 
non-trivial SPT ground state, for any partition of the chain, 
acquires an n-fold degeneracy. Remarkably, we shall verify 
that these unitary mappings, via Eq. O endow the parent 
Hamiltonians of non-trivial SPT chains with the required n- 
fold degeneracy of the energy spectrum per edge. Once the 
unitary transformation Eq. is derived, we can obtain the 
most general form of the SPT ground state Eq. 0 by al¬ 
lowing (3 ^ 0 without changing either the gapped nature of 
the many-body energy spectrum or the topological properties 
of the ground state. That this is true can be seen perturba- 
tively: moving away from the /3 = 0 limit with (3 « 1 and 
Eg{s) a local function in Eq. (2.151, amounts to adding small 
local symmetry-preserving perturbations to the gapped (3 = 0 
theory, which therefore, cannot immediately destroy the SPT 
phase. Based on this we expect that, for 0 < /? < (3 C , the 
ground state Eq. <o describes a class of many-body SPT 
ground states adiabatically connected to the /? = 0 limit. 


and the SPT state 


^ZaXZa ) = W Z 2 xZ 2 I $Z 2 xZ 2 ) 

= vLe eiWz2xZ2is) I s ) 


(3.5) 


Le* 1 


2tVs/2 


Ei @j (1 8 3 S j + 1)/2| 


»>■ 


The unitary transformation Eq. endows the state 

Eq. (3.51 with an amplitude e lSj for every domain wall be¬ 
tween neighbor spins in the many-body configuration | s ). We 


seek to find 9j for which Eq. (3.5 i describes the SPT ground 
state. 

Due to the product state nature of Eq. ( |3.2[ ) and the pairwise 
entanglement induced by the unitary mapping Eq. ( |3.4| >, one 
effortlessly finds that, for any partition E*, the reduced density 
operator obtained by tracing over one of the subsystems is 
given by the 2x2 matrix 


A. Z 2 x Z 2 SPT states in d = 1 

We consider a spin chain, with an even number N s of sites, 
divided into even and odd sublattices and Z2 x Z2 symmetry 
generated by 


c(l) _ 
- 


n 

j6even 


Xi 


o(2) _ 

— 


n 

odd 


X; 


(3.1) 


1 


Pi = 


1 cos (0i) 
cos ( 0i ) 1 


(3.6) 


For 0i = 0 the above density matrix has a single non-zero 
Schmidt eigenvalue. The existence of 2 degenerate Schmidt 
eigenvalues (equal to 1/2) is verified for 


di = ±tt/ 2. 


(3.7) 


In the following, { X , Y, Z} denote the three Pauli matrices. 

We start by examining a trivial spin chain in the special 
limit where it is represented by a product state. 


AC 


$7/ 


>=n 


j =1 


, 1 

>i 2 n */ 2 


Eld 


(3.2) 


Moreover, imposing that the unitary transformation 
Eq. (|3.4| > commutes with either of the Z 2 symmetries in 
Eq. (|3.1|> yields the final form 


= I I e *f 


n 


'• f (1-^2 j Z 2 j + l) 


(3.8) 


where Xj \ —1). = 


and { | s ) } represents the 2 


AC 


many-body spin states in the eigenbasis of Z operators: 


’ > = 


ij ), for Sj ± 1. We recognize Eq. (3.21 as 


the ground state Eq. ( |1.2| ) in the 8 = 0 limit. It has the parent 
Hamiltonian 


H; 


Z2 XZ2 


N s 


— u o/Ai~Xj) 

3-1. 


(3.3) 


which is an RK Hamiltonian of the form Eq. (2.21 with (3 = 0 
and u! StS i = x'o > 0, where the sum in Eq. (|2.2a 1 extends over 
pairs of states (s, s') which differ by a single spin flip. 

The product state form of the ground state Eq. ( |3.2[ > implies 
that its Schmidt decomposition, with respect to any partition 
Ej between sites i and 1 + I of the lattice, contains a single 
Schmidt eigenvalue (equal to 1). 

Now, for 0j £ R, we introduce the unitary mapping 


W! 


Z2 XZ2 


n 


= lie 


(1 Zj Zj+ 1)/2 


(3.4) 


One then finds 


x z 2 x 2 j w Z2XZ2 - X 2j Z 2 j —1 Z 2j+1 

r—1 


X, 


21+1 W Z 2 XZ 2 - X 2 j+1 z 2j Z 2j+2 . 


(3.9) 


The operator Wz 2X z 2 Xj can be regarded as a 

modified Pauli spin operator (since it is obtained from Xj by a 
unitary transformation) which is “dressed” by the domain wall 
operator Z :) _ 1 Z J+ \ with support on the opposite sublattice. 

So, under the unitary transformation Eq. (3.81, the SPT 
Hamiltonian at zero correlation length is 


%Z 2 XZ 2 — 


2 X ^ 

N s 


H- 


Z 2 xz 2 W Z2XZ2 


— w 0 


E ( X ~ Z 3-l X 3 Z i+1 


(3.10) 


l=i 


We note that the model Eq. ( 3. 10[ ) has been constructed in 
Ref. (21 using the concept of decorated domain walls, while we 
have arrived on it by appealing to the notion of entanglement 
spectrum via the unitary mapping Eq. 0- 
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Ground state degeneracy can be easily attested by studying 
this model with open boundary conditions, where 


N,. 


/T/open _ 




i=2 


(3.11) 


The fact that the above Hamiltonian commutes with Z\ and 
Zj\r s implies that there are 2-fold degenerate states associated 
with the left and the right edges corresponding to states with 
si = ±1 and SjVs = ±1. Thus the universal properties of 
the SPT state studied here are encoded in the unitary mapping 
Eq. ((378). 

With the unitary transformation Eq. ( |3.8| ), we can obtain 
the more general form of the SPT ground state Eq. ED b y 
allowing (3 ^ 0 without changing either the gapped nature of 
the many-body energy spectrum or the topological properties 
of the ground state. In this regard, the Hamiltonian Eq. ( |3.10| > 
is a particular exam ple of a larger class of SPT models de¬ 


scribed in Eq. (2.151, with the phase factors e 




2 («) 


given 

by acting with the unitary transformation Eq. ( |3.8| >. on the 
state | s). 


which is an RK Hamiltonians of the form Eq. (2.21 with (3 = 
0, u) SiS ' = wo > 0, where the sum in Eq. (|2.2ai extends over 
pairs of states (s, s') which differ by a single Z 3 spin flip. 

We now consider, for 0 :i £ R, the unitary transformation 


w z 3 xz 3 =n ex p 


and the SPT state 


■n Y' 1 ( a j a j+iy 


(w a — 1) (w a — 1) 


(3.16) 


I ^Z 3 xZ 3 ) = Wz 3xZ3 I $Z 3 xZ 3 ) 

s 5K«E e, “ w * l i s > 


3 N ./2 

1 

3^/2 


E ex p 


» 1 _ u} a(s 1 + 1 s j ) 

^ 1 J (w a - 1) (w a - 1) 

7 a —1 v ' v ' 


|s). 

(3.17) 


B. Z3 x Z3 SPT states in d = 1 

We consider a spin chain, with an even number N s of sites 
and Z 3 x Z 3 symmetry generated by 


$£= n n, 


n 


j Geven j G odd 

where, at each site j , we consider the operators 

/0 0 1^ 

Tj = 1 0 0 

\0 1 0; 


(3.12) 



(3.13) 


satisfying r? = <j| = 1, r| <jj Tj = u> a j, where w = e l27r / 3 . 

Let | Si ), for Si = 0,1,2, be th e eigenstates of ai: 
<7i | Si) = w Si | Si ). As in Sec. (Ill Ai, we start with a triv¬ 
ial Z 3 x Z 3 paramagnet described by the ground state 


I $7/ 


>-n 


) + | lj ) + | 2j ) 

V3 


1 


3^/2 


E I s )’ 


14 


(3.14) 


where { | s ) } represents the 3 ■ many-body spin states in the 


eigenbasis of a operators. We recognize Eq. (3.141 as the RK 


ground state Eq. (1.21 in the (3 = 0 limit. It has the parent 
Hamiltonian 


N s 

#Z 3 xZ 3 = Wo E [ 2 - ( T J + T j ) 


3 =1 


(3.15) 


T he un itary transformation Eq. ( |3.16[ ) endows the state 
Eq. (3.171 with an amplitude e*for every pair of neighbor 


spins j and j + 1 for which Sj 7^ .s J+ i (mod 3) in the many- 
body configuration | s). 


Due to the product state nature of Eq. (3.141 and the pair¬ 


wise entanglement induced by the unitary mapping Eq. (3.161, 
one effortlessly verifies that, for any partition Sj, the reduced 
density operator obtained by tracing over one of the subsys¬ 
tems is given by the 3x3 matrix 


Pi = 


1 m 
m) 1 
M) m 



(3.18) 


where f(8) = |(1 + 2 cos ( 9 )). For 0j = 0 the above density 
matrix has a single non-zero Schmidt eigenvalue. The exis¬ 
tence of 3 degenerate Schmidt eigenvalues (equal to 1/3) is 
then verified for 


P= 1,2. 


(3.19) 


Moreover, imposing that the unitary transformation 
Eq. ( 3. 16| > commutes with either of the Z 3 symmetries in 
Eq. p.121 yields the final form 
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w; 


(p) 

Z3 XZ3 


n . 2 np ^ 

eX p * 


1 - ( a h-i °2jY 


“ (iU a — 1) (uj a — 1). 

a=l v 7 v 7 


exp 


. 27175 1 - (o'lj a 2j+i) 


3 ^ (uj a - 1) (u“ - 1) 

a—1 v 7 v 7 


P=l,2. 


(3.20) 


Eq. ( 3.20| > establishes two unitary mappings between the triv¬ 
ial SPT ground state and the two non-trivial SPT ground 
states. 


Moreover we find that 


w: 


(?) 


r 2 j (Wg } xz3 ) — T 2 j (&2j-l 0’2 :J - + 1 ) P , 


-1 


wgU^+r (W^ XZ3 )' =T 2j+1 {a\ 3 a 23+2 y, 


(3.21) 


so that the operator 77 gets “dressed” by a domain wall like 
operator - carrying Z 3 charge p - 1,2- with support on the 
opposite sublattice. 

The SPT Hamiltonians in the (3 = 0 limit then read 


y(p) — W 
rt z 3 xz 3 — w : 

= w 0 


(p) 


-1 


#z 3 xz 3 (wg ) xZ3 ) 

J2{ 1 ~ T2 i ( CT 27-1 4j+lf 


+ H.c.j + W 0 ^2 { 1 - T 2 j +1 (fj'ij a 2 j +2 ) p + H.c.j , 


(3.22) 


forp = 1, 2. 


As in Sec. ( Ill A| ). degeneracy of the ground state energy 
manifold can be checked by placing this system in an open 
chain, in which case, one finds that the SPT Hamiltonians 
commutes with 07 and thus implying a 3-fold degenerate 
state associated to the left and right edges. Thus the univer¬ 
sal properties of the SPT state studied here are encoded in the 
unitary mapping Eq. ( |3.20| i. 

With the unitary transformation Eq. ( 3.20| > we can obtain 
the more general form of the SPT ground state Eq. {□} by 
allowing (3 7^ 0 without changing either the gapped nature of 
the many-body energy spectrum or the topological properties 


of the ground state. In this regard, the Hamiltonian Eq. ( |3.22[ ) 
provides a par ticula r example of a larger class of models de¬ 


scribed in Eq. (2.151, with the phase factors e 


*Wz, 




given 

by acting with the unitary transformation Eq. ( |3.20[ ). on the 
state Is). 


x Z ?1 SPT states in d = 1 


We generalize the findings of Secs. (III A) and (IIIBi to 


Z„ x Z„ SPT states (see the Appendix for definitions and 
useful formulas). We arrive at n — 1 Z„ x Z„ symmetric 
unitary transformations, labeled by p 1 ,n — 1, 


W. 


(p) 


x Z„ 


n . 27t p y^ 

exp *- > 

n 


71 — 1 


1 - a 2j) a - 

(w a - 1) (w a - 1). 


exp 


27 Tp \—-v 

- « - > 

n ' 


71 — 1 


1 - (alj <J2j+i) a ■ 

(w a - 1) (d) a - 1). ’ 


(3.23) 


It can be shown, using Eqs. (A31, (|A4|> and (A5 1, that 


w: 


(p) 

7j n x Z ? 


T 2j x z „) = t 2 j (a 2 j—i 4j+iV 


w: 


zjxz„ r 2j + l 




— T 2j+ 1 (a 2j <J 2 jj r2 ) p , 

(3.24) 


forp = 1 , ...,n — 1 . 

The Z n x Z„ SPT Hamiltonians in the /3 = 0 limit are of 
the same form as Eq. (|3.22 1 with p = 1,..., n — 1. We notice 


r 


that this class of one dimensional SPT models has been stud¬ 
ied in Ref. [26] without reference to the connection between 
unitary transformations and ground state entanglement spec¬ 
trum that we are exploring here. Moreover, in our formalism, 
we do not need to be restricted to the /? = 0 limit for, with 


the unitary transformation Eq. (3.23!, we can construct a large 


class of SPT-RK Hamiltonians of the form Eq. (2.151 with 
the phase factors e 1 WZn xZ ™ given by acting with the unitary 
transformation Eq. (3.231 on the basis state | s ). 
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IV. SPT STATE IN ONE DIMENSION WITH 
TIME-REVERSAL SYMMETRY 

We now study the one dimensional bosonic SPT states pro¬ 
tected by time-reversal symmetry which, being an anti-unitary 
symmetry, requires conditions Eq. ( 1 .8b[ ) to be satisfied. 

The action of time-reversal symmetry shall be represented 
by the anti-unitary operator 0, 


ee _1 = i, © 2 = i, 

where we work with the representation 


(4.1a) 


due to the fact that, with periodic boundary conditions, 
(1 — Sj Sj+ 1)/2 is an even number. Moreover, comparing 
Eqs. ( |4.5[ ) and ( |4.7[ ) with the results obtained in Sec (III A i, we 
conclude that, for any partition E,; between sites i and i + 1, 
there are 2 degenerate (equal to 1/2) Schmidt eigenvalues and 
that the SPT parent Hamiltonian for the ground state Eq. ( |4.4[ > 
reads 


n b 


Ujt = W Z T Hjt w z / = w 0 ^2 (1 + Zj- 1 Xj Z j+ 1 ) . 

3 =i 

(4.8) 


0 - ( n X J ) K ’ (4.1b) 

3 

with K denoting the complex conjugation operator. 

A one dimensional chain of N s sites in the trivial Z/ insu¬ 
lating phase can be described by the Hamiltonian 

N s 

(1-Xj) , (4.2) 

3=1 

which has the time-reversal symmetric ground state 

1%) = yQi 51 l s >> (4.3) 

S 

where { | s ) } represents the 2 Ns many-body spin states in the 
eigenbasis of Z operators: Zj \ Sj ) = Sj s 7 ), for s 7 = ±1. 
We recognize Eq. ( |4.3[> a s the ground state Eq. ( | 1 ,2[ ) in the 
(3 = 0 limit and Eq. ( |4.2[ > as an RK Hamiltonian of the form 
Eq. \22\ with (3 = 0 and u) a y = wo > 0, where the sum in 
Eq. ( |2.2a[ > extends over pairs of states (s, s') which differ by a 
single spin flip. 

We now obtain the unitary mapping W z t, which, applied 
on the trivial ground state Eq. ( |4.3[ ), yields the non-trivial SPT 
state 

I % > = W z t I % ) = —^ X] e iW ^ {s) | a), (4.4) 

S 

wherein the action of W z t on the many-body basis {| s )} is 
defined as 


^N 7j T\s) = e W ^ s) \s). (4.5) 


Invariance of the SPT ground state | 'f z r ) under time- 
reversal symmetry Eq. ( |4.1[ ), according to Eq. ( |1.8b| >, implies 

W Z T (s) = -WjT ( s) + 2t rx integer. (4.6) 


In light of the discussion in Sec. ( pH} , if we assume for 
W z t an ansatz that leads to the entanglement of nearest neigh¬ 
bor spins, than the only transformation (up to a trivial global 
phase) consistent with Eq. ( |4.6| > has 

W T?M = \ E 0 --8383+1) (4.7) 

3 


Since this is essential ly the same model as the Z2 x Z2 spin 
chain discussed in Sec (III Ai, we conclude that the Z^ sym¬ 
metric SPT Hamiltonian has 2-fold degeneracy of the ground 
state per edge. 

With the unitary transformation Eqs. ( |4.5[ > and ( |4.7| i we can 
obtain the most general form of the SPT ground state Eq. 
by allowing (3^0 without changing either the gapped na¬ 
ture of the many-body energy spectrum or the topological 
properties of the ground state. In this regard, the Hamilto¬ 
nian Eq. (|4.8|» provides a particular example of a larger class 
of models described in E q. ( |2. 15| >, with the phase factors 

e Wz % ^ S> given in Eq. (4.71. 


V. SPT STATES IN TWO DIMENSIONS 
A. Z2 SPT states in d = 2 

The classification of bosonic SPT states in two dimensions, 
protected by Z„ symmetry, asserts the existence of n gapped 
phases of matterP®From the point of view of the one dimen¬ 
sional edge, these n phases can be distinguished by the way 
right- and left-moving degrees of freedom transform under the 
symmetry: while in the trivial phase right- and left-moving 
modes carry the same Z n charges, their Z„ charges are differ¬ 
ent for each of the other n — 1 non-trivial SPT states.® 15 As 
a consequence, in each of these n — 1 non-trivial SPT states 
the edge states cannot be gapped and symmetry preserving at 
the same time. 

Levin and Gu have constructed a model on a triangular lat¬ 
tice that describes a Z2 SPT paramagnet in two dimensions. 20 
The Z 2 symmetry is generated by 

Sz 2 = Y[X 0 , (5.1) 

3 

where {X, Y,Z} denote the three Pauli matrices and 

Z 3 I S 3 ) = S 3 I S 3 )’ f0r S 3 = ±L 

In Ref. [20] the two kinds of paramagnetic ground states are 
represented by 

l $z =) = ^72 Z I s )’ ( 5 - 2 ) 

S 

for the trivial paramagnet, and 


(5.3) 
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for the non-trivial paramagnet. N s denotes the number of 
sites of the triangular lattice and L(s) £ N counts the number 
of loops that separate domain-wall regions for each one of 
the 2 Ns many-body spin configurations | s) (see Figs, [l] [ 2 ] 
an d 0 . Thus while the coefficients in the expansion Eq. ( p.2[ 
all have the same sign, in the non-trivial SPT ground state 
Eq. £3}, the sign of the coefficients depends on the number 
of loops L(s). This non-trivial (real) phase factor was shown 
in Ref.|20]to be responsible for the universal properties of the 
Z 2 non-trivial paramagnetic state. 

From Eqs. ( [5T2] > and ([53]), we identify the ground states of 
Ref. [20] with Eqs. ( | 1 .2[ > and ( | 1.1 1 ) in the limit [3 = 0, and the 
unitary mapping Wz 2 between the two classes of paramag- 
nets, defined by its action on the basis {| s )}, to be 


ical properties of the ground state give rise to the edge physics 
studied in Ref. J 6 

t t t t 



W Z2 | s) =e lW ^ s) \s) = e i7ri(s) | s). (5.4) 


Given that the unitary transformation defined in Eq. ( |5.4| > 
encodes the universal properties of the two dimensional Z 2 
SPT paramagnet, we can move away from the “infinite tem¬ 
perature” limit (3 = 0 studied in Ref. (f2Ql) by allowing a 
larger class of Z 2 SPT paramagnets with ground state given by 
Eq .(O and microscopic Hamiltonians given by Eq. (|2.15|l. 


B. Z 2 x Z 2 SPT states in <1 = 2 


t t t t 



FIG. 1. (Color online) A particular many-body configuration | ) 

of spins Z (red), where the loops are defined along the domain 
walls and have support on the dual lattice. As in Ref. [20J in order 
to determine the projection of the Z 2 symmetry on the edge, we take 
the outer spins in a reference state where they are all pointing up. 


In Ref. M the classification of two dimensional bosonic 
SPT states with Z n x Z m symmetry was studied from the point 
of view of the effective one dimensional edge degrees of free¬ 
dom. The non-trivial SPT edge states were there shown to be 
associated with non-trivial 3-cocycles in the group cohomol¬ 
ogy. However, no microscopic theory for the two-dimensional 
bulk states has been investigated. In this section, we construct 
a microscopic theory for the Z 2 x Z 2 bosonic SPT states using 
the formalism of RK Hamiltonians and show how the topolog- 


FIG. 2. (Color online) A particular many-body configuration | ) 

of spins (blue). 


t t t t 



FIG. 3. (Color online) A particular many-body configuration 
| ) of the spins Zj 12 - 1 = Z j 1 ' Zj 2> (green). Note how the di¬ 

rection of the arrows at every site is consistent with Figs.[I]an d0 


In addressing SPT states with Z 2 x Z 2 symmetry, we allow, 
at each site j of the lattice, two spin species represented by 
Pauli operators { X f 5 , y} x) , zf ] } and { xf ] , y} 2) , zf }, 
where the Z 2 x Z 2 symmetry is generated by 

cl 1 ) _ TT vl 1 ) e( 2 ) __ TT y ( 2 ) 

“II X :i ’ “ll A i ■ (5.5) 

j j 

We immediately realize that there should exist at least four 
kinds of paramagnets, which are obtained by stacking two de¬ 
coupled Z 2 symmetric systems. These four states of matter 
are parametrized by two integer numbers, pi,p 2 £ { 0 , 1 }, 
which contribute to the expansion of the ground state with 
the (real) phase factor exp{i7r [pi L(s^) + p 2 L(s {2) )]}, 
where L{s^) and L(s^) count, respectively, the number of 
loops defined by domain wall configurations formed by spins 
of species 1 and 2, as in Figs.[l]an 00 

In order to obtain Z 2 x Z 2 symmetric states that go be¬ 
yond the tensor product of two Z 2 symmetric states, we real¬ 
ize that, at each site j , we can consider an independent Ising 
spin Zj 12 ' 1 = Z ’j 1 -* zj 2 ' with eigenvalues sj- 12 " 1 = sj 2 \ 

We then propose 
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2d 

Z2 XZ2 - 


^ gZTT [pi i(s ( 1 ) )+p 2 i(s < 2 ) )+Pl2 L(S (12) )] 

^ S W s( 2 ) 


(5.6) 


as a description of the 8 classes of Z 2 x Z 2 symmetric ground 
states parametrized by 3 binary indices pi,P 2 ,Pi 2 £ {0,1}. 
We adopt the notation | ) < 8 > | s*- 2 -*) = | sC 2 ' 1 ). When 

Pi = P2 = P12 = 0, Eq. ( |5, 6 [ > describes the trivial ground 
state 

I^2XZ 2 ) = ^ l s(1) ; s(2) >> (5.7) 

s(!) a( 2 ) 


which is the tensor product of two trivial Z 2 paramagnets. 

Eq. ( |5.6| > corresponds to the choice of ground state in the 
“high temperature” limit /3 = 0 and the SPT ground states 
can be obtained from the trivial one, Eq. ( |5.7| >, via the uni¬ 
tary transformation W| d xZ2 , whose action on the many-body 
basis { | sW; s^ ) } reads 


W 2d 

w Z 2 xZ 2 


s (l) ;s (2) ^ = e zW| 2 d XZ2 ( 


s (t). s (2) ^ _ gZTT [pi Z(s (1, )+P2 Z(s (2) )+P12 Z(s (12) )] 


S W:sM 


)■ 


(5.8) 


In Eqs. ( |5.6| l and ( |5 . 8 [ >. in order to simplify the notation, we 
omit the indices p\, p 2 and p i 2 on which the SPT ground states 
and the unitary transformation depend upon. 


1. Parent Hamiltonians of the Z 2 x Z 2 state in d = 2 


describes a trivial paramagnetic system and has Eq. as 
its ground state. It is the sum of two RK Hamiltonians of the 
form Eq. (2.2 1 with /? = 0, u SiS i = loq > 0, where the sum in 
Eq. (2.2a I extends over pairs of states (s, s') which differ by a 
single spin flip. 


The Hamiltonian 




M 


1 - xf ) 


(5.9) 


n 


2d 

Z2 XZ2 


= w: 


2d 


Hr 


2d 


r2d 


XZ2 Z2 XZ2 \ Z2 XZ2 / 


7 2d 


Then applying the unitary transformation W z xZ 
Eq. ( |5.8| ) to the Hamiltonian Eq. ( |5.9[ ) yields the microscopic 
model realizing the 8 classes of Z 2 x Z 2 SPT paramagnets: 


= Wo ^ |i - X-i) e i7r ( p 1 +p i2)+ i 5 [pi M" 4 1 ) ) + p i 2 ( 1 -^ 12> 4 12) )]} 

3 

+ Wo ^ |l_ J^(2) e i7r (P2+Pl2)+i 5 J2 <ee f ;j> [p 2 (1 -z‘ 2) Z f <2) )+Pi2 (l-z' 12) x'! 2) )] | 


(5.10) 


where the notation < ££': j > denotes that the summation 
in the exponent runs over pairs of nearest neighbors £ and 
£', which themselves belong to the 6 nearest sites around a 
given site j of the triangular lattice. When p-\ 2 = 0 and 
Pi = P2 = T Eq. ( 5.1 ()[ > describes, up to an overall energy 
shift, two decoupled Levin-Gu Hamiltonians.® On the other 
hand, the case pi 2 = 1 encodes strong entanglement between 
the two Ising systems. 


2. Edge symmetry transformations of the Z 2 x Z 2 state in d = 2 

The physical properties of the edge of two dimensional 
bosonic SPT states can be understood from the chiral ac¬ 
tion of the symmetry on the edge model ® 5 E 1 S |33 which. 


on a microscopic scale, originates from a non-onsite symme¬ 
try realization on the edge degrees of freedom . 9 20 For the 
Z 2 x Z 2 SPT states considered here, the non-onsite symme¬ 
try transformations on the edge were shown in Ref.[l 6 ]to be a 
manifestation of the non-trivial type-II cocycles of group co¬ 
homology. 

The ground state Eq. is a simple product state and, as 
such, the edge states are disentangled from the bulk. On the 
other hand, in the other SPT ground states whose expansion 
coefficients depend upon the domain wall configurations, the 
effect of flipping an edge spin can change the number of loops 
L , implying an entanglement between edge and bulk degrees 
of freedom. Therefore by studying the properties of the edge 
states, we can gain information about the nature of the bulk 
states and vice versa. 
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We have found, working directly with the bulk SPT wave 


functions Eq. ( |5.6[>, that the projection of the symmetry trans¬ 
formations Eq. (|5.5|> onto the boundary spins is 


ci 1 ) 

‘- > Z 2 ,edge 


n* 


(1) in (P1+P12)+* f Z'?, , 1 )+p 1 2(l-^ (12) 


(5.11a) 


3(2) 

‘- , Z 2 ,edge 


n* 


(2) i n (p 2 +pi 2 )+i f [p 2 (l —• Z j 2) -Zj + i) +P 12 (l-zj 3 * * * * * * * * 12) zj +1 )] 


(5.11b) 


where the product in Eq. ( |5.1 1 1 > is taken over the spins on the 
edge with periodic boundary conditions. Remarkably these 
are the same non-onsite symmetry transformations studied in 
Ref. 021 which shows that Eq. ( |5.6[ > gives the correct repre¬ 
sentation of the ground state of the two dimensional bosonic 
SPT states protected by Z 2 x Z 2 symmetry. As is clear from 
Eq. ( |5 .11 1 >, the non-onsite form acquired by the edge sym¬ 
metry in the non-trivial SPT phase implies that domain wall 
configurations on the edge carry projective representation of 
the Z 2 symmetry, i.e., upon flipping any given spin configu¬ 
ration twice, there is a factor of (—l) p for every domain wall. 
Since the total number of domain walls is even, it follows that 
Eq. ( |5.11[ ) is a faithful representation of the Z 2 symmetry on 
the edge. 


3. A refinement of the group cohomology classification in 

terms of cobordisms has been recently proposed to de¬ 
scribe three dimensional bosonic SPT states with time- 

reversal symmetry.® It could be worthwhile to study 

these time-reversal symmetric phases of matter using 

the formalism of RK Hamiltonians, in particular, to un¬ 

derstand if there is any special meaning attached to the 

global unitary transformation that connects the trivial 
and the non-trivial SPT ground states, which could shed 

new light on the reason behind the inability of the group 

cohomology classification to account for these phases 
of matter. 
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els whose ground states encode the universal properties of 
bosonic symmetry-protected topological states. Although we 
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tonians for bosonic SPT states in two dimensions with 
discrete symmetry Z n , for n > 3. In contrast to the 
Z 2 case, the expansion of the paramagnetic state in the 
ordered basis contains more than one type of domain 
walls, suggesting that the coefficients in the expansion 
of the SPT ground state in the | s ) basis may depend 
on orientable loops. It remains an open question how 
these domain wall fluctuations can coherently give rise 
to SPT ground states. 

2. We expect that this formalism may offer a path to 
constructing microscopic models for other classes of 
bosonic SPT states in three dimensions, as well as in¬ 
teresting examples of bosonic SPT states in two dimen¬ 
sions protected by Z 2 x Z 2 x Z 2 symmetry, which are 
classified by non-trivial Type-III cocycles!-® 
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Appendix A: Z„ operators 


At each site j of a lattice we consider operators ay and r, 
represented bynxn matrices 


/1 0 0 0 \ 

0 u> 0 0 


\0 0 0 w ""7 


(Al) 


/0 0 0 ... 0 1 \ 

1 0 0 ... 0 0 

0 1 0 ... 0 0 

0 0 1 ... 0 0 

\: 0 0 ... 1 0 / 


satisfying 


= r; = 1 . 


J 3 ‘3 


J 3 5 


(A2) 


(A3) 
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where lu = e 1 and lu = e %2 n . 

The following identities hold for p = 0, n: 

J 2 np 

exp < i - 

l n 

(A4) 


n- 1 ^ [ a J a j' 

2 + ^ w“ - 1 

a= 1 






exp 



a— 1 
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